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TRANSPORT-COLLAPSE SCHEME FOR SCALAR 
CONSERVATION LAWS INITIAL AND BOUNDARY VALUE 

PROBLEMS 

DARKO MITROVIC AND ANDREJ NOVAK 


Abstract. We extend Brenier’s transport collapse scheme on the Cauchy 
problem for heterogeneous scalar conservation laws and initial-boundary value 
problem for homogeneous scalar conservation laws. It is based on averaging out 
the solution to the corresponding kinetic equation, and it necessarily converges 
toward the entropy admissible solution. In the case of initial-boundary value 
problems, such a procedure is used to construct a numerical scheme which 
leads to a new solution concept for the initial-boundary value problems for 
scalar conservation laws. We also provide numerical examples. 


1. Introduction 

The subject of the paper is the construction of new numerical method for Cauchy 
and initial-boundary problems for scalar conservation laws. The method is a gen¬ 
eralization of the transport-collapse scheme introduced in [5]. A consequence of 
the analysis of the scheme is a new solution concept of the initial-boundary value 
problem for scalar conservation laws, and this is the most important contribution 
of the paper. 

In order to introduce it, let fl C R"^ be a bounded smooth domain and R+ = 
[0,oo). We consider 

dtu + diVx/(M) = 0, (t, x) G R”*" X O, (I) 

w|t=o = wo(x), (2) 

mIr+xOH = (3) 

where / G C^(R;R‘^). If not stated otherwise, we assume that uq G L^(R‘*) C 
RV(R‘^), ub G L[q^(R+ X dfl) n i3V(R+ x dfl). We also assume that 

a < uq, us < 6 for some constants a < b. (4) 

A typical problem described by (1), (2), (3) arises e.g. in traffic flow models. 
Namely, if we aim to describe a flow on a finite highway (required to model on and 
off ramps) we need to use boundary conditions [22]. For instance, optimization of 
travel time and cost between two points can be obtained by controlling incoming 
and outgoing car densities [2]. 

Nevertheless, it is clear that the boundary conditions cannot be prescribed if the 
characteristics corresponding to equation (7) and emerging from the boundary leave 
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n. This means that one needs to introduce a new concept defining what conditions 
the unknown function u should satisfy in order to be a solution to (1), (2), (3). 
This was first done in [4] where the existence of strong traces at the boundary of 
solutions is assumed; see [1, 20, 23]. The weak formulation which does not require 
existence of strong traces was later proposed by F.Otto [19] and the corresponding 
numerical method was developed in [24]. The concept is further extended in [21] 
in a more general setting (on manifolds necessarily implying that the flux depends 
on x). Let us recall it here. 

Definition 1. A function u G is said to be the weak entropy solution to 

( 1 ), ( 2 ), (3) if there exists a constant L G R such that for every fc € R and every 
non-negative (p G Cc(R+; R^), R+ = R"^ x R^^, it holds 


[ {\u-k\+dtip+ sgn_^_{u-k){f{u) - f{k))V^ip) dxdt (5) 

-I- / \uo — k\+ip{0,-)dx +L / (p\uB — k\+d^{x)dt > 0, 

Jr'^ jR+xan 


and 


/ {\u - k\-dtp> + sgii_{u - k){f {u) - f{k))V^(p)dxdt ( 6 ) 

+ / \uo — l\-(p{0, ■)dx + L / ip\uB — k\-d'y{Tc)dt > 0, 

Jr<^ JR+xdn 

where 7 is the measure on dO,. 

As it comes to the refinement of the latter concept, we actually base it on an 
interesting observation from [ 21 ] roughly stating that if the characteristics emerging 
from {t = 0} X n hit the boundary then the corresponding boundary value should 
not affect the solution. Thus, we are going to construct the definition of solution 
so that it involves somehow only those parts of the boundary which essentially 
influence on solutions. 

To be more precise, assume that we are dealing with the flux depending on x 
i.e. / = /(x. A). Denote by 

5“ = {x e 9D : {fx{x,X),v)<0 a.e. A G/}, 

where I contains all essential values of the functions ub and uq (i.e. of appropriate 
entropy solution u), and v is the outer unit normal on dVL. The set S~ actually 
consists of all points such that all possible characteristics from that point enter 
into the (interior of the) set D. Therefore, for every x G S~ , the trace of the 
corresponding entropy solution is actually ub(x). 

Similarly, for 


5'“*'= {x G i9D : (/)^(x, A), z7) > 0 a.e. A G /}, 
all possible characteristics issuing from x G 5+ leave the set D, and u_b(x) does not 
influence on the weak entropy solution u to (1), (2), (3). 

However, both sets S~ and S'+ can be empty since for some A G /, it can be 
(/)^(x, A), < 0 and for other A G / we could have (/)^(x. A), i7) > 0. Therefore, 

in order to refine former arguments, we need to rewrite considered conservation 
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laws SO that we can more accurately take into account behaviour of the flux / with 
respect to A. A natural choice is the kinetic formulation to (1) since it includes 
the variable A in a desired way. Before we introduce it, let us recall the Kruzhkov 
entropy admissibility conditions for (general, heterogeneous) scalar conservation 
laws. 


Definition 2. A bounded function u is called an entropy admissible solution to 


+ diVx/(t,x,M) = 0 (7) 

with the initial conditions (2) if for every convex function V € (7^(R), every A G R 
and every non-negative (p € it holds 

JJ [V{u)dt(p+J f'^(t,x,v)V'{v)dv ■Vip+ J diY^f{t,x.,v)V'\v)dv(p\dxdt 


R+xR'* 


( 8 ) 


+ [ K(uo(x))</5(0,x)dx < 0; 


Remark 3. In the case of the heterogeneous equation (7), we shall consider the 
Cauchy problem. As we shall see, the heterogeneity causes significant technical 
challenges. It is possible to overcome them in the case of the boundary problem as 
well, but we believe that this would be unnecessary complication which would hide 
main ideas of the new initial-boundary concept. In the numerical examples at the 
end of the paper, we shall take a space dependent flux in order to show how the 
method works in general situations. 


Equivalent and more usual definition of admissible solution is given by the 
Kruzhkov entropies K(A) = |n — A|, A G R, and it states that a bounded func¬ 
tion u is called an entropy admissible solution to (7), (2) if for every A G R it 
holds 


9t|u- A| -bdivx[sgn('u- A)(/(t,x,u) - /(t,x. A))] -f sgn(M - A)divx/(t, x. A) < 0 

( 9 ) 


in the sense of distributions on I?'(R'^), and it holds esslimt^o |u(t, x)—uo(x)|dx = 
0. Roughly speaking, by finding derivative with respect to A in (9) one reaches to 
the kinetic formulation provided below (see e.g. [11, 14, 16] for different variants). 


Theorem 4. [8] The function u G C([0, oo); nL^^((0, oo); L°^(R‘^)) is the 

entropy admissible solution to (7), (2) if and only if there exists a non-negative 
Radon measure m(t,x. A) such that x R'^"*'^) < oo for all T > 0 and such 


{ 1, 0 < A < M 

— 1, It < A < 0 , represents the distributional solution 

0 , else 

to 


dtx + div(x,A) [F{t, X, A)x] = dxm{t, x. A), {t, x) G R+ x R"*, 

X{\u{t = 0,x)) = x(A,uo(x)), 

d 


where E =(/),,- X] d^jfj)- 


( 10 ) 

( 11 ) 
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In the next section, we shall provide properties of the function x- 
Remark that through the kinetic concept, one reduces the nonlinear equation 
(7) on the linear (so called kinetic) equation (see Theorem 4). However, derivative 
of a measure hgures in the equation (see the right-hand side of (10)) and it has 
one more variable which is usually called kinetic or velocity variable. Due to the 
former reason, the kinetic equation is not convenient for numerical implementation. 
Nevertheless, if we neglect the derivative of the measure, and then average out the 
solution to the obtained linear equation with respect to the kinetic variable, we 
obtain entropy solution to the considered problem. Such a procedure is proposed 
in [5] for Cauchy problems corresponding to equation (1). One of the aims of the 
paper is to extend the transport-collapse scheme [5] for the initial value problem 
for heterogeneous scalar conservation laws. 

The power of the method to be presented is in its ability to transform nonlinear 
problem into linear. Linear scalar conservation laws are easy to solve numerically 
since there are a lot of robust numerical schemes available. The cost of that ’’trans¬ 
formation” in practical computing is adding one more dimension (see (10)). 

Moreover, we shall use the transport-collapse techniques to construct the bounded 
function u satisfying the following definition. 

Definition 5. We say that the function u G L-(R+ X 17; [a, 6]) is a weak entropy 
admissible solution to (1), (2), (3) if for every k G R and every non-negative 
C'i(R$;R+) it holds 



kl+dt(p + sgn+(u - k){f{u) - f{k)) Vdxdt + [ 


|iio - k\+(p{0, ■)dyi 


( 12 ) 


pb—a p 

J J ip\X — k — a\ + {f'{X — a),iy{x))x{^,UB{t,x) — a)d'y{x)dtdX > 0, 


R+xao 

{f'{\-a),P(x))<0 


and 



(|u — k\-dt(fi + sgn_(u 


k){f{u) - f{k))\7:^(f) diidt+ [ \uo - k\-(fi{0,-)dx 

jR'i 


(13) 


(p\X — k + b\-{f'{X + 6), v(x))x(A, UB^t, x) — b)d'y{x)dtdX > 0. 


1—b J 

R+xan 

</'(A-l-h),i7(x))<0 


It is not difficult to see that if u satisfies conditions of Definition 5 then u also 
satisfies Definition 1. This will be proved in the last section. 

Let us briefly comment Definition 5. The first two terms on the left-hand sides 
of (12) and (13) are standard in the entropy admissibility concept (compare with 
Definition 1 and (9)) and they are related to the behaviour of the solution u in the 
interior of D and on t = 0. The final terms on the left-hand sides of (12) and (13) 
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simply say that when the characteristics enter Q (i.e. when the angle between the 
normal z? and /'(A) is greater than tt/2, i.e. when (/'(A_i7(x)) < 0) then we shall 
take the boundary data into account. Remark that we shifted a solution by a in 
(12) and by b in (13) since then u — a > 0 and u — & < 0, respectively, implying 
x(A,u(t,x) — a) > 0 and x(A,u(t,x) — 5) < 0. This enabled us precise control of 
the behaviour of the solution at the boundary (see the last section). 

Finally, let us remark that work in the field of numerical methods for conser¬ 
vation laws is rather intensive. Most of the papers deal with Cauchy problems 
for conservation laws (scalar conservation laws or systems; see e.g. classical books 
[10, 15] and references therein). As for ( 1 ), ( 2 ), (3), there are not so many results 
since the interest for this kind of problem has arisen relatively recently. We mention 
[3, 24] and references therein. For results in the case of systems, one can consult 
[18] where one can also hnd thorough overview of the subject. 

The paper is organized as follows. In Section 2, we shall prove convergence of 
the transport-collapse scheme for initial value problems corresponding to (7). In 
Section 3, we shall introduce a transport-collapse type operator for (1), (2), (3), 
and the proof of its convergence toward the entropy solution. 


2. Transport collapse scheme for the Cauchy problem for 

HETEROGENEOUS SCALAR CONSERVATION LAW 

The transport-collapse scheme is based on tracking of characteristics of conser¬ 
vation law (10). In the homogeneous case (i.e. when the flux is independent of 
(<,x)), the characteristic have very simple form x — and it is significantly 

easier to analyse them than in the case when the flux is (t, x)-dependent. Thus, 
this section represents a non-trivial generalization of the method from [5]. 

Let us first introduce assumptions on the flux / from (7). We assume / = 
(/i) ■ • • j /d) € C'^(R(J.’^^) and for some b > a 

f(t,x,a) = /(<,x,6) = 0, 

and a < uq < b for the initial condition uq. Latter conditions provide the maximum 
principle for the entropy admissible solution to (7), (2). More precisely, the entropy 
admissible solution u will stay bounded between a and b. 

Let us now state properties of the function x- 

Proposition 6. [5, page 1018] It holds 

a) Vu, V G L^(R'^ ) such that u> V x(A, u) > x(A, v); 

b) Vu G L^(R‘^), V(? G ^““(R^^ X R), it holds 
II X(A, u) 5 '(x, A)dxdA = / (/“ ^(x, X)dX) dx; 

In particular, if g = G'^ and G{a) = 0, then ff x(X, u)g(x, X)dxdX = 
I G(x,u)dx 

c) TV{u) = I TV{x{\-))dX. 

The idea of the transport collapse scheme for the initial value problem (7), (2) 
is to solve problem (10), (11) when we omit the right-hand side in (10): 


dth V div^xiP{t, X, A)/i] = 0, h\t=o = x(A, uo(x)). 


(14) 
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The solution of this equation is obtained via the method of characteristics. They 
are given by 

I ^ = fx, x|t=o = xo, 


^ = - T, A|t=o=Ao. 


I, f=i 

For later purpose, we rewrite this system in the integral form 


(15) 


X =yio + J* 

A = Ao - /o X) fj {t\ X, X)dt'. 
i=i 


(16) 


The solution to (14) has the form 

h{t,x,X) = x(Ao(t,x, A),Mo(xo(t,x, A))). (17) 

To avoid proliferation of symbols, denote for tp = {(pi,... ,tpfi) : —>■ the 

norms 


II ^xV^II oo 


sup sup 

Ax>0 


|(pfc(x+ Ax) - (/9fe(x)| 
|Ax| 


||‘/?||oo = sup sup |<Pfc(x)|, 
l<k<d xSi'R.'^ 


(18) 


where | • | denotes the Euclidean norm. We have the following properties of the 
characteristics. 


Proposition 7. The characteristics Xq = Xo(t, x, A) and Aq = Ao(t,x, A) satisfy 
the following continuity properties: 

l^xl := |xo(t, X + Ax, A) - Xo(t, x, A)| (19) 

< |Ax| (^1 + y maxIlVx/U^'r, A)||oodty ■ 

I^^aI := |Ao(t, X + Ax, A) - Ao(t, x, A)| (20) 

< |Ax| f max||Vxdivx/(t',-, A)||oodt', 

Jo ^ 

where the norms are given by (18). 

Proof: From (16), we have 


x = Xo(t,x,A)+ [ ,x,X)dt' 

Jo 

X + Ax = xo(t,X + Ax,A) + / fxit',x + Ax,X)dt'. 

Jo 


By subtracting those equations, we obtain: 

|xo(t,x +Ax, A) -xo(t,x, A)| 


<|Ax|+ f max||/(^(t',x +Ax, A) -/j^(t',x, A)||oodi' 

Jo ^ 

< |Ax| + |Ax| [ maxIlVx/K^',-, A)||oodi'- 
Jo ^ 


(21) 
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This proves (19). Inequality (20) is proved analogously. It holds 

rt d 


A = Ao(t,x +Ax, A) - [ da:-fj (t',X + Ax, X)dt', 

Jo 

A = Ao(t,x,A)- [ 

Jo 


and it is enough to subtract the last two equalities, and to follow the procedure 
from (21). □ 

Let us now define the transport-collapse operator T. 

Definition 8. The transport collapse operator T{t) is defined for every u G 

by 

T{t)u{yL)= j x(Ao(t,x, A),M(xo(t,x, A)))dA. (22) 


ft satisfies the following properties which are the same as the ones from [5, 
Proposition 1], 

Proposition 9. It holds for every u,v G L^(R‘^) 

a) u < V a. e. implies T(t)u < T(t)v a.e; 

b) / r(t)it(x)dx = f u(x)dx; 

c) the operator T(t) is non-expansive 

\\T{t)u-T{t)v\\L^^d) < ||M-u||il(R,d), 

and, in particular, ||T(t)M||Li(Rd) < |L||Li(R<i); 

d) TV{T(t)u) < (1 -I- Cit)TV{u) tC 2 , where TV is the total variation and 
Cl and C 2 are appropriate constants depending on the -bounds of the 
flux f; 

e) |jT(t)rt — u|| j;^i(Rd) < C 2 TV{u)t -\-tCi for the constants Ci and C 2 from the 
previous item; 

Proof: Item a) directly follows from the definition of the transport collapse op¬ 
erator T{t). 

As for the item b), for every fixed t > 0, denote by Z = (t,x, A) characteristics 
from (15). Notice that, since div(x,A)^ = 0, it holds 


dZ{t,XQ,Xo) 

9(xo,Ao) 


= 1 . 


Therefore, according to Proposition 6, 


(23) 


T{t)u{x)dx = 


/R-i+l 


x(Ao(t,x, A),M(xo(t,x, X))dxdX 


(24) 


/xo(t,x,A) = yA 
VAo(t,x, A) = r]) 


/Rd+l 


Xiv,u{y)) 


det 


dZ{t,xo,Xo) 


d(xo,Xo) 


dydr] = J u{y)dy. 


Item c) now follows from a) and b) according to the Crandall-Tartar lemma 
about non-expansive order preserving mappings [7, Proposition 3.1]. 
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Let US now prove item d). We have 
\T{t)u{x + Ax) — T(t)u(x)|dx 




= y I X+ Ax, A), u(xo(t, x+Ax, A)))-x(Ao(xo(t, x. A), u{xo{t, x, A)))dA|dx 

R'i R 

< J k(^o(L x+Ax, A), u(xo(t,x+Ax, A)))-x(Ao(xo(t, x. A), u(xo(t, x, A)))|dx(iA 


Rrf+i 


We next write Xo(t, x + Ax, A) = Xo(<, x. A) + 7^x(L x. A) and Ao(t, x + Ax, A) = 
Ao(t, X, A) + TZ\{t, X, A), where TZ^ and TZx are estimated in (19), and introduce the 
change of variables Xo(t, x. A) = y, Xo{t, x. A) = 77 (keep in mind (23)). We obtain 

I |T(t)M(x + Ax) — T(t)u(x)|c?x 

< I \x{v + T^\,u{y + Tly,))-x{v,u{y))\dydr] 
jR'i+l 


< 


[ Ixid + 'Jlx,u{y + TZ^j) - x{v,u{y + 'Jl^))\dydr] 
jR'i + l 

Ixid^uiy + n^)) - x{v,u{y))\dydr] 


jR'i+i 

< ||7^A||ooTl/(x) + l|7^x 


TV{x{vM-)))dv = 4||7^A||oo + \\TZ^\UTV{u), 


'R 


since the characteristics are of C^-class, TV (x) = 4, and since Proposition 6, item c) 
holds. Remark that in the case when u>0 we actually have x(A, u) = sgn_|_(M — A) 
and in that case TV{x) = 1. Having in mind Proposition 7, we conclude the proof 
of d). We remark that 

Cl = 4inax||Vx/(t, •, A)||oo, C '2 = inax || Vxdivx/(t, •, A)||oo- 

It remains to prove item e). Using (16), as in to the proof of item d), we have 

||T(f)u-w||ii(Rd) < [ |x(Ao(t,x, A),'u(a;o(t,x, A))) - x(A,M(x))|dxdA 
Jr<‘+i 

= [ lx(A + 7^A,^i(x + 7^x))-x(A,^i(x))|dAdx 
jR‘i+1 

< CitTViu) + C2t, 

which immediately gives e). □ 

We also need the following result. 

Proposition 10. For any smooth positive test function (p, any u G L^(R) such 
that a < u < b, and convex Lipschitz function U : R —)■ R, we have 

J {V{T{t)u) — V{u))(x.)(p{x.)dx < J J Bv{t',x,u{x))Wipdxdt' (25) 


m divx/(t, X, X)V''{X)dXdt' + oft), t —)• 0 

. 
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where Bv{t,K.,u) = X)V'{\)d\, and o{t) depends only on the L°°-bound 

ofu. 

Proof: Remark first that for any fixed (t,x), from the definition of the function 
X, it follows for any -function G 


/ G'(A)x(Ao(t, X, A), uiMt, X, A)))dA = ^(-l)'=G(wfe) - G(0), (26) 

k=0 

where the increasing sequence (uJk), k = 0,... ,2p, belongs to the set Nul = {A e 
[a, 5] : Ao(t,x, A) = M(xo(t, x, A))} (since the entropy solution to (7), (2) takes 
values in the interval (a, 6)). Indeed, for almost every (t,x) G x R, it holds 
Nul = {A S [a, 6] : Ao(t,x, A) = u(xo(t, x, A))} = {ojo, ■ ■ ■ ,co 2 p}. In the intervals 
(wfcjWfe+i) and (wfc+ijWfc) the function A i—)■ Ao(t,x, A) — u(xo(t, x, A)) has different 
signs, and we can assume 

Ao(t,x,A) > u(xo(t,x,A)), A € (uj 2 s,co 2 s-i-i), s = 0,...,p; 

Ao(t,x,A) < u(xo(t,x,A)), A G (w 2 m-si,W 2 m), m = 0,...,p. 

To be more concise, recall that we assumed a = 0 and 6 > 0. According to 
definition of the kinetic function y, we see that x(Ao(t, x. A), xo(t, x, A)) = 1 for A G 
(a; 2 m+i, W 2 m), m = 0,...,p, and x(Ao(t, x. A), Xo(t, x. A)) = 0 for A G (a; 2 s, W 2 s-i-i), 
s = 0,... ,p. From here, (26) immediately follows. 

Remark that the set has odd cardinality since the multivalued solution is ob¬ 
tained by continuous transformation from the graph of initial value [5, page 1016]. 
Moreover, due to the mean value theorem, the following relation holds for any 
convex function V (see e.g. [9, p. 51]): 


2p 2p 


From (26) and (27), it follows 


(27) 


V(T(t)u(x)) = V{ x(Ao(i,x, A),M(t,xo(t,x, A)))dA) (28) 

< ^(-l)'=R(a;fc) / V'(A)x(Ao(t,x,A),u(xo(t,x,A)))dA-hV(0). 

k^O 

Moreover, it holds 


F(r(x)) = I V(A)x(A, u{x))dX - ViO). (29) 

Subtracting (43) from (28), we reach to 

V(T{t)u{x)) -V{u{x)) < j y'(A) (x(Ao(t,x, A),'u(xo(t,x, A))) - x(A, u(x))) dA. 

(30) 
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We have from here 

J (y iT{t)u{x)) - V(u(x))) <p(x)dx (31) 

< jj (^'(A)x(Ao(t,x, X),u{xo{t,x,X)))-V'{X)x{X,u{x)))^{x)dxdX 
= JJ V'iXoit, X, A))x(Ao(l, X, A), u(xo(t, x, A)))((/j(x)-vj(xo(<, x, X)))dxdX (32) 
+ JJ {V'{X)-V'{Xo{t,x,X)))x{Xoit,x,X),u{xo{t,x,X)))(p{x)dxdX (33) 

+ ( J J V' {Xo{t,x, X))x{Xo{,t,x, X),u{xQ{t,x, X))) (p{xQ{t,x, X))dxdX 

- JJ V'{X)x{X,u{x))ip{x)dxdXy (34) 

The two terms from (34) cancel according to (23). Indeed, using the change of 
variables from (23), with the notation from (24), we conclude 


JJ V'{Xo{t, X, A))x(Ao(f, X, A), u(xo(t, x. A))) i^(xo(t, x, X))dxdX 
= JJ V'{r])x{v, u{y))>p{y)dydr]. 


Let us now consider the term from (32). Using the Taylor formula 


JJ V'iXoit, X, A))x(Ao(t, X, A), u(xo(t, x, A)))(v3(xo(t, x, A)) - ipix))dxdX (35) 

= JJ V'iXoit, X, A))x(Ao(t,x, A),m(xo(<,x, A)))(xo(t, x. A) - x) • V(^(xo(t,x, X))dxdX 

+ ^ JJ V'iXoit, X, A))x(Ao(t,x, A),M(xo(t,x, X)))D'^ipix) (xo(t,x, A) -x)'^dxdX, 

where x is a point belonging to a neighbourhood of xoit, x. A). To further estimate 
the latter term, we expand the function f'xit', x. A) into the Taylor expansion around 
Xq and take (16) into account: 


x-xo(t,x,A)= f f'xit' ,x,X)dt' = f f'xit', xoit,x,X),X)dt' (36) 

Jo Jo 

+ (xo(t,x. A) - x) • / X/a;f'xit',x,X)dt'= f f'xit',xoit,x,X),X)dt'+ Oit'^), 

Jo Jo 

since clearly Xo(<, x. A) — x = Oit) and /J f'xit', x, X)dt' = Oit). Inserting this into 
(35) and applying the change of variables from (23), we conclude using item b) 
from (6): 

JJ V'iXoit, X, A))x(A, uix))iipixoit, X, A)) - (pix))dxdX (37) 

= f Bvit',x,uix))W(pdxdt'+ Oit^). 

Jo 

To deal with the remaining term from (33), we shall expand the function V' into 
the Taylor series around Aq. We have 
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- ^'(^o))x(Ao(i,x, A),M(xo(t,x, A)))(p(x)dxdA (38) 

= jj V"{Xoit, X, A))(A - Ao((t, X, A)))x(Ao(t, x, A), m(xo(<, x, A))) (p(x)rfxdA 
+ C>(||A - Ao((t,x, 

Applying the procedure as in (36) and having in mind (16), we reach to the estimate 


Ao(t,x, A) — A = — 



(39) 


If we notice that IIA - Ao((t, x, A))||L(supp(vp)x(a.b)) = ^om (38) and (39), we 
conclude 


- ^'(Ao))A(Ao(t,x, A),M(xo(t,x, A)))(p(xo(t,x, A))dxdA (40) 

pt p pu d 

= -J J J ^/j(t',xo(t,x,A),A)v5(x)dAdxdt' + C>(t^) 


Combining (31), (37), and (40), we conclude the theorem. □ 

A consequence of Proposition 9 and Proposition 10 is the following theorem: 

Theorem 11. Denote 

S„(t)u = (1 - a)T{-fu + aT{-f+^u, (41) 

n n 

where 

fceN, as [0,1). (42) 

n 

For each initial value Uq G L^(R‘^) such that a < Uq < b, the unique entropy 
solution of (7), (2) at time t is given by the formula 

u(t,-) = — lim Sn{t)uo. 

n—>-oo 

Proof: First, fix an arbitrary t > 0. Consider the sequence of functions Unit, •) = 
Snit)u. We aim to prove that the sequence (M„(t, •)) is strongly precompact in 
L^(R'^). To this end, we shall use the Kolmogorov criterion stating that a functional 
sequence bounded in L^(R‘^) is strongly precompact in L^(R‘^) if it is uniformly 
L^(R‘^) continuous. In other words, we need to prove that 

a) \\unit, •)IIli(R‘^) ^ C for every n G N and some constant (7; 

b) for any relatively compact K CC R'^, any e > 0, there exists Aa; > 0 such 
that ||M„(t,x +Ax) - it„(t,x)|jii(ad) < e. 

Item a) follows from Proposition 9, item c) (we take u = 0 there). 

As for the item b), we shall use (recursively) property d) from Proposition 9. Tak¬ 
ing into account definition of the total variation and form of the sequence (M„(t, •)), 
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simple calculations show that (with the notations from Proposition 9) 

TV{un) < {l + Ci-)TV{T{-T-\uo)) + C2- < ... 

n n n 

- (l ^ TV{u) + ^ ^ (^1 + j < exp{Cit)TV{uo) + 2C2{t + t^). 

Since uq G -BV(R‘^), this immediately implies L^-equicontinuity of the sequence 
(un(t, •)). This means that for every fixed t > 0, we can choose a strongly converging 
subsequence (not relabelled) (it„(t, •)) of the sequence (un(t, •))■ By taking a dense 
countable subset i? C R+, we can choose the same converging subsequence (it„(t, •)) 
for every t G E. 

Now, by the continuity property given in item e) from Proposition 9, we conclude 
that the subsequence (un{t, ■)) strongly converges in C([0, Tl; L^(R‘^) for every T G 
R+ toward a function u G C([0, T]; ^^(R'^). 

Now, we need to check that u satisfies the entropy admissibility conditions. First, 
notice that for every t, as n —>■ oo, it holds that a —>■ 0. Thus, it is enough to notice 
that the main part of the transport-collapse operator given by T{-)^u -G u as 
n ^ oo along the previously chosen subsequence and to consider 



{V{T{-y+^u) - F(T(-)Ju))(^(x)dx 
n n 


(43) 


(25) /-(t + l)*/" f f 

< / / Bv{t',x.,T{ — yu{x))Vipdxdt' + 0(t/n). 

^ Jjt/u JR-i n 


Now, we simply let n —>■ oo and keep in mind arbitrariness of t to infer that the 
function u satisfies the entropy admissibility conditions from Definition 2, a). 

Remark also that this implies convergence of the entire sequence given by (41) 
due to uniqueness of entropy solutions to (7), (2). □ 


3. Boundary value problem 

In this section, we shall consider boundary value problem for homogeneous scalar 
conservation law (1) on the domain D, which is a bounded simply connected open 
smooth subset of R'^. 

In order to simplify the presentation, we shall assume that a = 0 in (4), i.e. that 
the solution to the considered problem is non-negative. In particular, this implies 
that the kinetic function y corresponding to such a solution satisfies 

x(A,u) = sgn+(u-A) > 0. (44) 

First, notice that the kinetic formulation from Theorem 2 still holds in the in¬ 
terior of R+ X However, we cannot use the method of characteristics from the 
previous section directly since the characteristics entering the boundary determine 
the value at the boundary. Nevertheless, since we are re-iterating the procedure 
after a short period of time (see (41)), we can modify the transport collapse scheme 
so that we take into account the boundary data. 
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Accordingly, recall that the kinetic reformulation for (1) has the form: 

5tx(A,u) +/'(A)divxx(A,M) = dxm+{t,x,X) (45) 

where m-^. is a non-negative measure. Assume that is an open set such that for 

some (T € (0,1), no two outer normals from dfl do not intersect in the set 

flo- = {x S : dist{x^ fl) < cr} i.e. in the set 

(i.e. we assume that has finite curvature). In order to augment (45) (with 
neglected right-hand side) with appropriate initial data, denote by i7(x), x € 
the unit outer normal on dfl passing trough the point x. We then extend the 
boundary data UB(t, x) for every fixed t > 0 along the normals p ( x ) in the set 
More precisely, we set for x S 11'^ = (slightly abusing the notation) 

ub(<, x) = ^^(t,Xg), for Xg S cin such that p(xo) = p(x). (46) 

Finally, introduce the function 

fO, X ^ 

= I u(t,x), X e (47) 

[u_B(t,x), X e V,cr\^ — , 

which is actually the extension of u along the normals v. If the function u does not 
depend on t, then we put u(x) instead of u(t,x), and wb ( 0 , x ) instead of UB{t,x) 
on the right-hand side of (47). Remark that we can rewrite the function .)(x) 
in the form 

w„(t..)(x) = u{t,x)Kn{x) + UB(t,x)Kn-(x), 
where ka is the characteristic function of the set A. 

Now, we are ready to introduce a modification of the transport collapse scheme 
from the previous section. Fix t > 0 and n S N. We neglect the right-hand side of 
(45) and, on the first step, we augment it with x(A, iCuo (x)) as the initial data. 

5th-I-/'(A)divx/i = 0, (48) 

h\t=o = (49) 

The solution to (48) is given by h(t,x, A) = x(A,Wuq(x— /'(A)t)) (since the char¬ 
acteristics of the equation have quite simple form; see [5]). We construct the ap¬ 
proximate solution Un to (1), (2), (3) by the following procedure: 

• 

Un{t',x) =T{t'/n){wuo{x)) := ( x(A, w„o(x -/'(A)t'))'^A, t'G(0,t/n]. (50) 

Jo 

• For fc = l,...,n — 1, we take 

Un{kt/n + t',x)= [ x(A,a;„„(fet/n,.)(x-/'(A)t'))rfA, t'G(0,t/n]. (51) 

Jo 

Remark that here, we have actually applied the transport collapse operator. Roughly 
speaking, the approximate solution in [0,t] x fl is given by the transport-collapse 
operator, while in [0,t] x the sequence (u„) is equal to the boundary data 
extended along the normals on dfl. 
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We shall show that the sequence (m„) strongly converges in L^{[0,t] x fl) along 
a subsequence toward a function u which represents the solution to (1), (2), (3) in 
the sense of Definition 5. In order to prove the later fact, we shall use the kinetic 
formulation similar to [21]. We introduce the following definition. 

Definition 12. We say that the non-negative function p_|_ G L°°(R+ x x R) is 
the kinetic super-solution to (1), (2), (3) if it satisfies the following equation for 
any ip G C,i(R5|^) and p G C'^(R) (recall that we assumed a = 0 in (4)): 


/r^xr 


p{k)p+{t, X, k) {dtp -f f'{k) VxP) dxdtdk 


(52) 


R Jo 


p{k)psgn^{\ - k){f{X),h!{x))xiX,UB{t,x))d'y{x)dtdXdk 


R+xan 
</'(A),!7(x))<0 


R JR+xQ 


p{t, x)p'{k)dm'^{t, X, k), 


for a non-negative measure m'^. 

We say that the non-positive function G L-(R+ X X R) is the kinetic 
sub-solution to (1), (2), (3) if it satisfies the following equation for any p G (7X11-+) 
and p & Cl (R): 


/R^ xR 


p{k)p- {t, X, k) {dtp + f'{k) Vxp) dxdtdk 
0 

RJ-b 


(53) 


[ [ [ p{k)psgii_^_{X-k + b){f'{X + b),ix{x))x 

JRJ-b J 

X x{^jUB{t,x) — b)d'y{x)dtdXdk 


R+xan 
(/'(A),j7(x))<0 


R JR+ xQ 


p{t,x)p'{k)dm {t,x,k), 


for a non-negative measure m~. 

The function p G L°°(R''' x x R) is the kinetic solution if it is kinetic super¬ 
solution and (1 — p) is kinetic sub-solution. 


The following theorem holds. 


Theorem 13. There exists the kinetic solution to (1), (2), (3) in the sense of 
Definition 12. 

Proof: We shall prove that the sequences of functions (sgn+(M„ — k)) (weakly) 
converge toward the kinetic super and sub solutions for the sequence (rt„) defined 
by (50) and (51). Since sgn_|_(it„ — fc) = 1 — sgn_(it„ — k) the kinetic super-solution 
will be the kinetic solution at the same time. 

First, remark that for every n G N, value of the function M„(t,x) for x G D is 
given by the transport collapse operator. 

Therefore, it is enough to consider behaviour of V{T{f)v) — V{v) for a convex 
function V whose special form will be chosen later, and for the function r; > 0 
playing the role of Un{ts, •) (recall that we have assumed that a = 0 implying that 
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our sequence of approximate solutions is non-negative) such that v(t, x) = usit, x), 
X e Accordingly, denote by Pn-/'(A)t(x) the unit vector to d{fl — f'{X)t), 

assume that we fixed t < a foi a given in (46), and consider for if G (fl^): 


(V(T(t)v)(x) - V(v)(x))f(x)dx 


(54) 


(30) 

< 


[ [ t^'(A)(x(A,u(x-/'(A)t))-x(A,i;(x)))(/?(x)(iAdx= (y = x-/'( 

[ [ V'{X)x{X,v{y))f{y + f{X)t)dydX- [ ( V'(X)x(X,v(x))f(x)dxdX 

Jo Jn-f'{\)t Jo Jn 

[ [ V'{X)x{X,UBiO,x))f{x + tf{X))dxdX 


(n-/'(A)t)\o 


r / 


t4'(A)x(A, v(x))ip{x)dxdX 


n\(n-/'(A)t) 


0 J {n-f'{\)t)nn 


V\X)x{X, v{x)){ip{x + tf{X)) - ip{x))dxdX. 


Now, since u > 0 it will also be x(A,z;) > 0 (see (44). Next, for a fixed fc G R, 
choose V{X) = 14+(A) = |A — k\+ in (54). After expanding the function ip in the 
Taylor expansion around x and taking into account that as t —>■ 0: 


i J g{x)dx^- j g{x){f{X),v{x))dx{x)>Q, 

(/'(A),^(x)><0 (/'(A).i7(x))<0 


we get (keep in mind that > 0) 


{V+{T{t)v){x) - V+{v){x)) ip{x)dx 


(55) 


< -t 


Vl{X){f{X),iJ{x))x{X,UB{t,x))ip{x)dxdX 


dQ 

(/'(A)T(x))<0 


+ t [ [ fWyiWx{Xv{t))'^{^)dxdX + o{t). 

Jn Jo 


Since for every n, the function Un(t,-) has the same properties as the function v 
from the above, we see that Un{t,-) satisfies (55). Therefore, as in the proof of 
Theorem 11 (more precisely relation (43)), we conclude that (m„) satisfies 
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k\+dt(p + sgn+(u„ 


k){f{un) - f{k)) Vxy>) dxdt 


(56) 


|mo - ^l+<p(0> ■)d.^ 

[ [ ip\X-k\+{f\X),d{x))x{X,UB{t,x))d'-f{x.)dtdX + O{-)>0. 

Jo J n 


/R'i 

rb 


R+x9n 


The left-hand side of the previous expression defines a non-positive distribution 
in (t, X, k) G R‘1 X R and thus, it is a non-positive measure. We denote it by 
m"(t,x,k). Having this in mind, we get after differentiating (56) with respect to 
k: 


I sgn^{un-k)+{dtip +f'{k)X/:^ip)dxdt+ sgn^{uo - k)ip{0, Jdx (57) 

R^ JR-^ 

J V5Sgn+(A - k){f{X),d{x))x{X,UB{t,x))dxix)dtdX 

R+xan 
(/'(A),j7(x))<0 

+ O(-) = 5fcTO+(t,x, fc), 

n 

where the above relation is understood in the sense of distributions in /c G R. 
Letting n —)■ oo in (57) along a subsequence for which both (sgn_|_(M„ — k)) and 
(m") weakly converge, we reach to (52). 

In order to get wanted relation for H(A) = H_(A) = |A — A:|_, remark that the 
function 



—b <w = u — b<0 


represents the weak solution to 

dtw + divx/(w -I- 6) = 0, 
with the initial and boundary data 


—b<wo=uo — b<0, —b< wb = ub — b < 0. 


If we apply the the transport-collapse procedure described in this section, then the 
corresponding sequence of approximate solutions has the form (wn) = {un — b) for 
{un) defined in (50) and (51). We can thus repeat the arguments from (54) to 
conclude (keep in mind that now x(A,'y(x)) < 0) 
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< 


{V_iT{t)iv)){x) - V_iv)(x))ip{x)dyi 

f 

J-b 


(58) 


V^^(A)x(A, ub(0,x) - b)ip{x + tf'{X + b))dxdX 


(n-/'(A+&)t)\n 


/° 

l-b 


VL (A)x(A, v{yi))tf{^)dyidX 


a\(n-/'(A+&)t) 


!-b J{n-f {x+b)t)nn 


^^(A)x(A, z;(x))((/j(x tf{X + b)) - ip{x))dxcdX. 


From here, as for V+, we obtain (53). □ 

Denote by p+ and p— the weak limits along appropriate subsequence of the 
sequences (sgn_|_(u„ — fc)) and (sgn_(M„ — fc)), respectively, defined in the proof of 
the previous theorem. It is not difficult to see that they satisfy conditions from 
[21, Definition 3.1]. This follows from the following theorem providing relation 
between Definition 5 and Definition 1 (from which the kinetic formulation given in 
[21, Definition 3.1] is derived). 

Theorem 14. A function u satisfying Definition 5 satisfies Definition 1. 


Proof: It is enough to notice 


|(r(A),p(x))|<L, 

where L is the constant such that |i/'||oo < L. Then, the last term in the left-hand 
side of (12) satisfies for every t > 0 (the other terms there are the same as the 
corresponding ones from (5)): 


nb—a n 

J J |A - fc - a|_^(/'(A - a),p(x))x(A,MB(t,x) - a)d 7 (x)dA| 


an 

(/'(A-a),j7(x))<0 


p po—a 

<L / / |A —/c — a|+x(A,MB(t,x) — a)dAd7(x) 

dan Jo 

= L |MB(t, x) — A:|c?x 
dan 


from where we conclude that u satisfies (5). The proof that u satisfies (6) is the 
same. □ 


Direct corollary of the previous theorem is existence and uniqueness of the function 
u satisfying Definition 5. 


Corollary 15. There exists the function u satisfying conditions of Definition 5 and 
it is unique. 


Proof: As we have already noticed, the functions p+ and p— constructed in 
Theorem 13 (see the comments after the proof of the theorem) satisfy conditions of 
[21, Definition 3.1]. Therefore, according to [21, Corollary 4.2], the function p_|_ has 
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the form p+{t, x, k) = sgn^(it(t, x) — k) and therefore, p-(t, x, k) = sgn._{u{t, x) — 
k) for some u G L°°(R+ x il). Now, it is a standard fare to conclude that u 
satisfies conditions of Definition 5 (it is the same as the proof of [21, Theorem 3.3]). 
According to Theorem 14 and the results from [19], we conclude that u is a unique 
solution to (1), (2), (3) in the sense of Definition 5. □ 

Corresponding numerical examples are given below. It is one-dimensional scalar 
conservation law defined on [0,0.5] x [—1,1] with the flux f{x,u) = i?e(x)(l — 
u){u -b 1 ) -I- 4:Hg{—x){l — u){u + 1 ), where is a standard regularization of the 
Heaviside function with e = 10“"^. In the first simulation boundary conditions are 
u\x=-i = 0, u\x=i = 1 and the initial condition is u\t=o — H^{x). In the second 
simulation boundary conditions are u\x=-i = 1, u\x=i = 0 and the initial condition 
is u|t=o = H^{-x). 




Figure 1. Cauchy problem (left) and boundary problem (right) 
with the initial condition uo(x) = Hg{x). 



Figure 2. Cauchy problem (left) and boundary problem (right) 
with the initial condition Uo(x) = Hg{—x). 
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